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Abstract. We use some of the recently released observational data to test the viability of two classes of minimally 
coupled scalar field models of quintessence with exponential potentials for which exact solutions of the Einstein 
equations are known. These models are very sturdy, depending on only one parameter - the Hubble constant. To 
compare predictions of our models with observations we concentrate on the following data: the power spectrum 
of the CMBR anisotropy as measured by WMAP, the publicly available data on type la supernovae, and the 
parameters of large scale structure determined by the 2-degree Field Galaxy Redshift Survey (2dFGRS). We use 
the WMAP data on the age of the universe and the Hubble constant to fix the free parameters in our models. We 
then show that the predictions of our models are consistent with the observed positions and relative heights of the 
first 3 peaks in the CMB power spectrum, with the energy density of dark energy as deduced from observations 
of distant type la supernovae, and bf with parameters of the large scale structure as determined by 2dFGRS, 
in particular with the average density of dark matter. Our models are also consistent with the results of the 
Sloan Digital Sky Survey (SDSS). Moreover, we investigate the evolution of matter density perturbations in our 
quintessential models, solve exactly the evolution equation for the density perturbations, and obtain an analytical 
expression for the growth index /. We verify that the approximate relation / ~ fijjj also holds in our models. 
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1. Introduction 

Recent observations of the type la supernovae and CMB anisotropy strongly indicate that the total matter- 
energy density of the universe is now dominated by some kind of dark energy or the cosmological constant A 
IjRiess fe al. 19981 IRiess 20001 IPerlmutter fe al. 19991 IRiess fe al. 2004 |l . The origin and nature of this dark energy 
remains unknown l|Zeldovich 19671 |Weinberg 19891 ICarroll 200l)l . 

In the last several years a new class of cosmological models has been proposed. In these models the standard cosmo- 
logical constant A-term is replaced by a dynamical, time-dependent component - quintessence or dark energy - that 
is added to baryons, cold dark matter (CDM), photons and neutrinos. The equation of state of the dark energy is 
given by wq = Pq/pq, Pq and pq being, respectively, the pressure and energy density, and —1 < wq < 0, which 
implies a negative contribution to the total pressure of the cosmic fluid. When wq = — 1, we recover a constant A- 
term. One of the possible physical realizations of quintessence is a cosmic scalar field, minimally coupled to the usual 
matter action l|Peebles fe Ratra 19881 ICaldwell fc al. 1 998) . Such a field induces dynamically a repulsive gravitational 
force, causing an accelerated expansion of the Universe, as recently discovered by observations of distant type la su- 
pernovae (SNIa) l|Perlmutter fc al. 19991 IRiess fc al. 19981 IRiess fc al. 2004 t and confirmed by WMAP observations 
( |Spergel fc al. 2003|). Accelerated ex pansion together with the strong observational evidence that the Universe is spa- 
tiallv flat ( Ide Bernardis fc al. 20001 |Spergel fc al. 2003| ) calls for an additional component and quintessence could be 
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responsible for the missing energy in a flat Universe with a subcritical matter density. Quintessence drives the cosmo- 
logical expansion at late times and also influences the growth of structure arising from gravitational instability. Dark 
energy could cluster gravitationally only on very large scales (> fOO Mpc), leaving an imprint on the microwave back- 
ground anisotropy IjCaldwell fc al. 19 98 ) ; on small scales, fluctuations in the dark energy are damped and they do not 
influence the evolution of perturbations in the pressureless matter ( IMa~S ^L_1999| . On the scales we are considering in 
the following, we assume that quintessence behaves as a smooth component, so that in our analysis the formation of clus- 
ters is due only to matter condensation, while the quintessence alters only the background cosmic evolution. The leading 
candidates for the dark energy as suggested by fundamental physics include vacuum energy, a rolling scalar field, and a 
network of slight topological defects l|Turner 2 000 Clin e 200T)l . Moreover, an eternally accelerating universe seems to be 
at odds with some formulations of the string theory l|Fischler fc al. 20 03 ). However, this is still controversial. Actually in 
the last few years it has been suggested that string theory could be compatible with the presently considered cosmologi- 
cal models (see for example, |Hellerman, Kalper fe Susskind 200T]lTownsend fe Wohlfarth 2003UGibbons fe Hull 200fl 
and references therein). This has stimulated a revival of interest in the exponential scalar field quintessence. In differ- 
ent scenarios, such exponential potentials can reproduce the present accelerated expansion of the universe and some 
predict future deceleration. Moreover, despite the criticism that exponential potentials require fine tuning, recently 
several authors have pointed out that the degree of fine tuning needed in this case is not greater than in other scenar- 
ios <|Cline 20011 IRubano fe Scudellaro 200T1 ICardenas fe al. 2002jl . 

In this work we show that the two models of quintessence for which general exact solutions of the Einstein equation 
are known are compatible with recent observational data, in particular with the power spectrum of the temperature 
anisotropy of the CMBR, the observations of type la supernovae, and the parameters of the large scale structure of 
matter distribution. As a first step we introduce a new parametrization of quintessence models considered by Rubano et 
al. and Rubano & Scudellaro (Ruba no fe al. 2 002 Ru bano fc Scudellaro 200T|l . which avoids the problem of branching 
of solutions highlighted in Cardenas & al., (2002). 

Since the physical features of these models will be extensively discussed in a forthcoming paper, in Sect. 2, we simply 
present the basic equations of the quintessence models used in this paper. In Sect. 3 the linear perturbation equation is 
solved for the two potentials and the growth of density perturbations is discussed; in Sect. 4 we show that predictions 
of our models are compatible with the observationally established power spectrum of CMB anisotropy, the SNIa data, 
and the results of estimates of the average mass density of the universe from galaxy redshift surveys. In Sect. 5 we 
discuss the possibility of constraining the equation of state of dark energy by age estimates of the universe. Sect. 6 is 
devoted to final conclusions. 



2. Model description 

In this paper we consider two quintessence models, with a single and a double exponential potential, for which exact 
analytic solutions are available. The discussion of the physical properties as well as the mathematical features of these 
models goes beyond the aims of this work; they are presented in Rubano & al. (2004). Here we will only give the 
basic relations. 



2.1. The single exponential potential 

We investigate spatially flat, homogeneous, and isotropic cosmological models filled with two non-interacting compo- 
nents: pressureless matter (dust) and a scalar field tp, minimally coupled with gravity. We first consider the potential 
introduced in Rubano & Scudellaro (2001), 



V(<p) cxexpj-y ' W 

For this potential the following substitution 

a = (u»)5 , (2) 

* = - (vf '<)) ■ (3) 

where a is the scale factor, makes it possible to integrate the Friedman equations exactly. Setting as usual a(0) = 
we have 

u = Uit, (4) 
v = v\t + v 2 t 3 , 
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Fig. 1. Time dependence of the scale factor a for the single exponential potential, for Ho = 0.97. 



where u±, v\ and vi are integration constants, so for a(t), we get 



a s {t) = u lV2 r{t 2 + — ) 

V2 



If by H we denote the time dependent Hubble parameter then 



4 

H =3t 



t 2 + ^- 



t 2 



V2 



(5) 



(6) 



To determine the integration constants ui, v\ and V2 we set the present time to — 1. This fixes the time-scale 
according to the (unknown) age of the universe. That is to say that we are using the age of the universe, to, as a unit 
of time. We then set ao = a(l) = 1, which is standard, and finally Ho — H(l). Because of our choice of time unit 
it turns out that our Ho is not the same as the Ho that appears in the standard FRW model. The two conditions 
specified above allow one to express all the basic cosmological parameters in terms of H - With these choices the whole 
history of the universe has been squeezed into the range of time [0, 1]. Moreover this model is uniquely parametrized 
by H only. Explicitly we have: 



t 



a a (t) = -[(3flo- 2)^ + 4-3^0], 
2 (2(3H - 2)t 2 + 4 - 3H ) 



H{t) 



% = 



3t {(3H - 2)t 2 + 4 - 3H ) ' 

(4 - 3H ) ((3ff - 2)t 2 + 4 - 3H ) 

[2(3if - 2)t 2 + 4 - 3H ] 2 
(3H - 2)t 2 (4(3Ho - 2)t 2 + 3(4 - 3H )) 



[2{3H - 2)t 2 + 4 - 3H } 2 
Therefore now the omega parameters of matter and the dark energy are 



^m = ^m(* = 1) 



n Vo = n v (t = i) 



2(4- 3 H ) 
- 9H 2 ' 

(3ff -2) (3g + 4) 
9 H 2 



The equation of state of dark energy evolves with time and the parameter w is given by 



1 



3(3^0-4) 



2 6(4 -3 H ) + 8 (3 H - 2) t 2 ' 
so that today we have 



(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

(13) 



w 



8-3H 
4 + 3H 



(14) 
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Fig. 2. Redshift dependence of the parameter w for the single exponential potential, for Ho = 0.97. 



In Fig. we show the time evolution of the scale factor while the redshift dependence of w is plotted in Fig. [21 
Asymptotically for t — ► oo, a(t) ~ t 4 / 3 and therefore in this model the universe is eternally accelerating and possesses 
a particle horizon. The relation between the dimensionless time t and the redshift z is given by 



2* 



[i 2 ((3H -2)< 2 +4-3ff )p 



(15) 



We shall compare the predictions of the above model with a flat cosmological model filled in with matter and the 
cosmological constant. In this case, the redshift dependent Hubble parameter is 



H{z) = H V«mo(1 + zf + (1 - Quo) , (16) 

where Ho is the standard Hubble constant. This class of models has two free parameters, Qmo and Hq, where 

So = 100 to kms-'Mpc" 1 . 

Let us assume that the age of the universe is 

to =7xlGy = 3.1510 16 7 s, 

where 7 is a constant to be determined by astronomical observations. With this definition it is possible to relate the 
value of Hq to the small h — -ffo/100 of the standard FRW model. It turns out that 



H = 0.1 h 7. 

If we accept that to = 13.7 Gy as given by the WMAP team ( |Spergel fc al. 2003| ) then H = 1.37h. 
2.2. The double exponential potential 

As the second example we consider the following double exponential potential: 



(17) 



V(<p) 



A exp 




(18) 



B exp 




2AB. 



where A and B are constants. As we see, if we take A, B > this model contains in a certain sense an intrinsic 
negative cosmological constant, but this does not lead to the re-collapsing typical in such cases. It turns out that this 
is an eternally expanding model, with alternate periods of accelerating and decelerating expansion. In the following 
we denote uj 2 — SAB. 

For this model exact solutions of the Einstein equations with scalar field exist, but their explicit mathematical form is 
rather complicated. General properties of these solutions are discussed in Rubano & al. (2004). Here we observe that 
if we use the same parametrization as in the previous case, e.g. taking the age of the Universe as a unit of time, the 
solution will depend on only two parameters: Hq and u). In order to get physically acceptable values of f2 v (i.e ft v > 0) 
we have to take u) < 2. However, in our analysis we use a more restrictive range, < uj < 1, as we will discuss in the 
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Fig. 3. Time dependence of the scale factor a for the double exponential potential, for Hq = 0.97, and lo = 0.1. 



next sections. We will see however that in this range of values of lo the matter density parameter Qmo is changing only 
slightly. In this case we obtain the following expressions for a, H and Q v , which are the main quantities we need: 



'(*) = 



H = - 

3 



(3H sin 2 lo - lo sin 2u)t 2 - (3H - 2) sin 2 wt 
2 sin 2 lo — lo sin 2lo 

( t (3H - 2) lo sin 2ut ^ 
2(3Hq sin to — lo sin 2lo) 



\ t2 (3H -2) 



(3^0 - 2) sin^ 



(19) 



(20) 



sin w — lo sin 



2cj/ 



- - 



(cot coso;t — sinu;£) 2 4- a; 2 sin 2 a;£(t 2 — 
(3H - 2) 



(3i? - 2) sin 2 
(3i/o sin 2 cj — w sin 2cj) 



(3i/o sin 2 lo — lo sin 2w) 



(3ifo-2Vsin 2u)t 



2(3i/o sin aj — 2cu cos a;) 



(21) 



where Hq and w are constants. Since we consider a flat model, CIm = 1 — Note that when i — > oo, fi^ — > const 
and in the generic case this constant is smaller than 1 and therefore, when t — > oo, Qm does not vanish, so at the late 
stages of evolution of this model dark energy and matter coexist. For large t, a(t) ~ t 2 ! 3 though O^, ^ and evolution 
of this model resembles the matter-dominated phase of the standard FRW universe, hence in this case the particle 
horizon does not appear. At the present epoch we get 



Q vo = %(t = l) = 

(3 H - 2) [3H (lo sin 2lo - 2 sin 2 lo) +2lo (sin 2lo - 2 lo)] 



9 H Q 2 (lo sin 2 lo - 2 sin 2 lo) 



n Mo =n M (t = i) 



' lo sin 2 lo + 3 Hq cos 2 lo — 2 lo 2 



3Hq ) (lo sin 2lo — 2 sin 2 lo) 
When to is small we obtain the following simple approximated formula for Q vo and O^o 
(3 H Q - 2) (30 (3 H + 4) - (9 H Q + 4) lo 2 ) 



27 Ho 2 (10 -lo 2 ) 



(22) 



(23) 



(24) 



2 (4 (30 - lo 2 ) - 3 i? (30 + lo 2 )) 
27 H 2 (10 -w 2 ) 



i2 M = o-7 rr2 7T7, ~ ■ ( 25 ) 



For lo in the range (0, 1), Eq. (I24|) shows that the dependence of f2 Vo on lo can be neglected. Actually, in analyzing 
the observational data we will set lo = 0.1 . In the double exponential potential model the relation between the 
dimensionless time t and the redshift z is given by 



6 



M.Demianski & al.: Two viable quintessence models. 




Fig. 4. Redshift dependence of the parameter w, for the double exponential potential, with Ho = 0.97 and u = 0.1. 



(lo sin 2uj — 2 sin 2 u>) 3 
1 + z = i '- 

[t 2 (u sin 2lu - 3H sin 2 w) + (3H - 2) sin 2 w t] 
which for small values of u> and t can be simplified to 

~ 1 2 
23t~3 

1 + z 



(26) 



(27) 



[A-3H Q -^(3H -2)u 2 ] 3 

Fig. |2| shows the time evolution of the scale factor a. As shown in Fig. 0] the equation of state for this double 
exponential potential evolves with time. The relation between Hq and the value of the Hubble constant in the FRW 
model is the same as in the single exponential case. 



3. Growth of density perturbations 

The equation describing evolution of the CDM density contrast, Sm = 8pm/Pm> f° r perturbations inside the horizon, 
is UPeebles 19801 IMa fc al. 19991) 



S M + 2H(t)S M - 4ttG P m6m = 0, 



(28) 



where the dot denotes the derivative with respect to time. In Eq. (|28|l the dark energy enters through its influence 
on the expansion rate H(t). We shall consider Eq. (|28() only in the matter dominated era, when the contribution of 
radiation is really negligible. 



3.1. The single exponential case 

For the model with the single exponential potential described by the Eq. QJ, the differential equation (|28|) reduces to 



d 2 6 M 4 [4-3ff -2(2-3ff )t 2 ] dS M 



2(A-3H )S 



M 



0. 



(29) 



dt 2 3f (4 - 3 H - (2 - 3 H ) t 2 ) dt 3t 2 [4 - 3H - (2 - 3 H )t 2 } 

Equation (|29|l is of Fuchsian type with 3 finite regular singular points, and a regular point at t — > oo; i.e. it is an 
hypergeometric equation, which has two linearly independent solutions, the growing mode 6+ and the decreasing mode 
6-. Solutions of this equation can be expressed in terms of the hypergeometric function of the second type 2^1 • We 
get 



5- oc - 



111 /3F - 2 
2' 3' 6' \3H -A 



and 

5 + oc ti 



2^1 



"1 


7 


11 


(3H 




3' 


6' 


6"' 


\3H Q 



(30) 



(31) 
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In the linear perturbation theory the peculiar velocity field v is determined by the density contrast IjPeebles 19801 
Padm anabha n 1993J) 



v(x) =H L f 6 M (y)-^^d 3 y, 

where the growth index / is defined as 
din 5m 



f = 



(32) 



(33) 



a is the scale factor. According to our conventions, using only the growing mode, we get 



In f -0.4 




-1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 

In n M 



Fig. 5. In / versus lnf^M in the single exponential potential model. 



4-3H + (3ff -2) t 2 
4 - 3 H + 2 (3 H - 2) i 2 
7 (3ff -2) t 2 

11 {3 H -4) 



| 2*1 


"1 


7 


11 


/3^o 












3 


6' 


y 


\3/r 





2F1 



4 13 17 fZH Q -2 



3' 6 ' 6 ' \3H - 4 



2*1 



(34) 



1 7 11 / 3g -2 > 2 
3' 6' 6 ' V3ff -4 



The growth index is usually approximated by / ~ SIm"- For ACDM models, a ~ 0.55 (see |Silveira &: Waga 1 994 
Wan g fc Steinhardt 1998|lLokas et al. 2004jl . In Fig. El we show the logarithm of / as a function of the logarithm of 
r^M for the single exponential potential. It turns out that a ~ 0.57 provides a good approximation to the model. In 
Fig. 6 we see that a can be considered a constant during the late stages of the universe evolution. 




Fig. 6. Redshift dependence of the exponent a in the single exponential potential model. We see that a can be considered 
constant during the cosmic expansion. 
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Fig. 7. Behaviour of the exact and approximate growth index / for the double exponential potential model. The solid line is 
the exact function and the dashed one the approximate function. 



3.2. The double exponential case 

For the model with double exponential potential described by Eq. I)18|l . the differential equation (|28J) becomes more 
complicated, it assumes the form 



dt 2 



2t (3Hq sin 2 to — lo sin 2lo) + lo sin(2o>i) (2 — 3 Hq) 



t 2 (3H sin 2 lo - lo sin 2 lo) + sin 2 (cot) (2 - 3H ) 



dS 



dt 



(35) 



4[3# 



sin 2 lo) + lo (sin 2lo — 2co)\ 



3 2 (3H - 2) sin 2 cot + 2t 2 (lo sin 2lo - 3H sin 2 



0. 



Eq. Ij35(l does not admit exact analytic solutions. However, since with our choice of normalization the whole history 
of the Universe is confined to the range t G [0,1], and since we choose co < 1, we can expand the trigonometric 
functions appearing in Eq. I|35[) in series around t = 1, obtaining an integrable differential equation, which is again a 
hypergeometric equation. For the growing mode we get 



<5+ oc f 3 2 F 1 



1 7 11 



(2 - 3H ) t 2 



3' 6' 6 ' 2[3ff -4-f(i?o-2V 2 ] 
We use the growing mode S + to construct the growth index /; according to Eq. I|33|l we obtain 



where 



A 1 



105 



(3flo - 4) 3 (2 - 3H )t 2 2 F 1 



x 33 



3#o-4--(# -2)o 



2F1 



4 16 17 (2-3# ) t 2 

3' Y' T' 2 [3iJ - 4 - §(fr - 2)w 2 ] 

1 7 11 (2 - 3ff ) t 2 

~3' 6' Y' 2 [3ff - 4 - §(iJ - 2)w 2 ] 



A 2 = i V 



(3F - 2) [w 2 (176 + 39(ff - 4)ff ) + (H Q - 2)(3H Q - 4)] 
3if l ■ 



(36) 



(37) 



(38) 



(39) 



--(4-3if ) [(3£f - 2)t 2 + 4 



In Fig. U| we see how accurate our approximate / is, in comparison with that obtained by numerical integration of 
Eq. (|35|l . In this case the growth index is approximated by / ~ !1m°, where a ~ 0.57. 



4. Observational data and predictions of our models 

The scalar field models of quintessence described in the previous sections depend on at most two arbitrary parameters 
that admit a simple physical interpretation. The cosmological model with single exponential potential of Eq. Q is 
uniquely parametrized by the Hubble constant Hq only, while the model with double exponential potential of Eq. I |18| l 
is parametrized by Hq and the frequency lo. We have shown, however, that the actual value of lo affects only slightly 
the most important quantities as, for instance, the density parameter fl Vo (see Eq. (1241) 1. Comparing predictions 
of this model with observations we will assume that co = 0.1. Using the WMAP data on the age of the universe 
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Fig. 8. Behaviour of the second derivative of the scale factor in the single exponential potential with Ho = 0.97, Q^g = 0.75. 
Please note that the transition from a decelerating to an accelerating expansion occurs close to z = 0.5, as predicted by recent 
observations of SNIa z t = 0.46 ± 0.13 JRiess fc al. 2004 I . 



to = 13.7 ± 0.2 Gy and on the Hubble constant H = 71 ± 5 kms^Mpc" 1 ( |Spergel fc al. 2003| ), from Eq.(17) we 
get that our H a = 0.97 ± 0.08. Once H and w are fixed our models are fully specified and do not contain any free 
parameters. In particular in the single exponential potential model we have tt Vo = 0.73 ± 0.09 and wq — —0.73 ± 0.06 
and in the double exponential potential model Vl Va = 0.74 ± 0.10 and u>o = —0.73 ± 0.07. We can now compare 
predictions of our models with available observational data to test their viability. In Fig. |H1 that with these values of 
Ho and h, the transition redshift from a decelerating to an accelerating phase in the evolution of the universe falls 
very close to z = 0.5, in agreement with recent results coming from the SNIa observations l|Riess fc al. 2004"^) . In the 
following we concentrate on three different kinds of observations: namely the measurements of the anisotropy of the 
CMBR by the WMAP team IjBennett fe al 20031 |Spergel fe al. 2003|), on the high z supernovae of type la, and the 
observations of the large scale structure by the 2dFGRS Team ( Hawkin s fc al. 2003|) . 



4.1. Constraints from CMBR anisotropy observations 

After the discovery of the cosmic microwave background radiation by Penzias & Wilson (1965) several experiments 
have been devoted to measuring the temperature fluctuations of CMBR. The recently released Wilkinson Microwave 
Anisotropy Probe (WMAP) data opened a new epoch in CMBR investigations, allowing strict tests of realistic cosmo- 
logical models. The WMAP data are powerful for cosmological investigations since the mission was carefully designed 
to limit systematic errors, which are actually very low (Ben nett fc al 2003"| . The WMAP measured the power spectrum 
of the CMBR temperature anisotropy, precisely determining positions and heights of the first two peaks. It turns out 
that the separation of the peaks depends on the amount of dark energy today, the amount at the last scattering, and 
some averaged equation of state. Assuming that the radiation propagates from the last scattering surface up to now 
in such a way that the positions of the peaks are not changed, the peaks appear at multipolc moments 

In = — -^-rm = uIa, (40) 

where r\ is the conformal time ( then r\ rec and 770 are respectively its value at the recombination and today), the 
speed of sound c s is assumed approximately constant during recombination, and Ia is the acoustic horizon scale. 
In models with quintessence Eq. (14011 should be modified and rewritten as (Doran & Lilley 2002, Doran & al. 2002 
IH11 fc Dodelson 20021 IBTi fc al. 200111 

l n = l A {n-(n), (41) 

where Cn = C + ^Cj here A£ is a general phase shift, whose form for the first three peaks is given in 
flDoran fc Lilley 2002| ), and C = ai(r rec ) a2 + f^ ec , where 

_ Pr( z rec) (42) 



■ ret / \ ' 

PM{Zrec) 
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The constants <xi and ai are fit parameters also furnished in Doran & Lilley (2002), Doran & al. (2002), and p r is 
the radiation density. It is possible to obtain an analytic formula for I a- Following Doran & Lilley and introducing 



( w o) = 



we get 



I A = 7TC S 1 



where 



yja rec {l - St™) 



F(Sl (po ,(w )) = 



dx 



o 



2o rec (l O^ ) 



n 



rll 



1 - 



-3<tuo) _j_ ^ro (1 - x) 



1 - fl 



(44) 



(45) 



(46) 



and a rec is the value of the scale factor at recombination. Using Eqs. (44), (|45|l and 1)41(1 in our single exponential 

potential model, we find that (w ) = —0.85, and we get the following values of I specifying the positions of the peaks 
l 



/i — 225 22 ? 

h 
h 



830tg . 



(47) 
(48) 
(49) 



These values are consistent with the observed positions of the peaks as measured by Boomerang 
(|de Bernardis fc al. 2000j) and WMAP flSpergel fc al. 2003| ): 



if 



BOOM 



iBOOM 



iBOOM 



= 539tf[, 



(50) 
(51) 
(52) 



WMAP 



22o.it2:| ; 



WMAP 



= 546 



+ 10 



(53) 
(54) 



In Fig. El we plot the CMB anisotropy power spectrum calculated for the single exponential potential model with 
H = 0.97, compared with the Boomerang data and WMAP l|de Bernardis fc al. 20001 |Spergel fc al. 2003| ). 
It is also interesting to calculate the relative height of the peaks. In particular Tii, which is the first peak amplitude 
relative to the COBE normalization ( that is, the height of the first peak is normalized with respect to the COBE 
result at I < 10), H.2, the relative height of the second peak with respect to the first, and H.3, the amplitude of the 
third peak with respect to the first. Using Camb l|Lewis et al. 2000}) with standard input parameters to evaluate the 
CMB power spectrum, we obtained Hi — 7.9, Ti.2 — 0.48, and H3 = 0.44, while a similar analysis performed on the 
Boomerang and Maxima data gives Hi = 7.6 ± 1.4, H 2 = 0.45 ± 0.04, and H 3 = 0.43 ± 0.07. 

As another complementary test we can use the WMAP result on the CMB shift parameter R — 1.710 ± 0.137, 
where ( |Hu fc Sugiyama' T996| ) 



R — Hq yj VLm q 



H(z) 



dz, 



(55) 



Zls 



1048[l + 0.00124(fVi 2 ) 



-0.738- 



[i + gi (n Mo h 2 y*], 



(56) 



where is the dimensionless baryon density and the functions g± and gi are given in Hu & Sugiyama (1996). According 
to Eqs. ||8J, (J3J, and l|15fl R can be written in terms of t as 



R = H oy /ShT (j:(t)\ to - F(t)\ t J), (57) 

1 We note that (wo) = —0.85 is fully compatible with the limit w < —0.78 usually accepted for the present equation of state 
of dark energy in quintessence models (Spergcl fc al. 2003| > . 
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Fig. 9. Behaviour of the theoretical CMB anisotropy power spectrum for the single exponential scalar field model with Ho 
compared with the data of Boomerang (empty boxes) and WMAP (filled boxes) 



0.97 



where T(t) is the dimensionless comoving distance: 



T{t) = 3(2t)3 



1 - 



(3ffo-2)t' 
3 Ho -4 



2^1 



1 1 

6' 3' 



(3H -2)t 2 
3 Ho -4 



H (f (3ff -2) + 4-3 J ff ) T 



(58) 



-,+0.07 
-0.07' 



so that R = 1.71^QQg, which is consistent with 



In our model, we obtain VL Mo h 2 = 0.13 ± 0.03, T = VL Mo h = 0.19 
the WMAP value. A similar analysis can be performed also for the model with the double exponential potential (1181) 
even if the calculations are more complicated and cannot be done by analytic procedures only. The first complication 
is naturally due to the presence of another parameter u>. However, as shown in Eg. ( 1241) . the final results depend only 
slightly on the value of to when picked from the acceptable range < w < 1. In analyzing the observational data we set 
u> = 0.1. Using Camb it is possible to calculate the locations of the CMB power spectrum peaks also for this double 
exponential potential model; we obtain 



h 

h 
h 



2241J 1 , 
825^. 



(59) 
(60) 
(61) 



In Fig. E3we plot the CMB power spectrum resulting from our model with the Boomerang and WMAP data. The 
evaluation of the relative height of the peaks gives identical results as in the single exponential potential case. In the 
same way we calculate (this time numerically) the observable quantity R according to the Eqs. i|55|) and Q56[l. and 
we obtain n Mo h 2 = 0.129 ± 0.05 ± 0.02, T = Q Mo h = 0.19ig;^lg;g|, so that R = 1.71±g:g|, which is consistent with 
the WMAP result R = 1.70 ± 0.07. We note that the second term in the errors takes into account the effect of the 
indeterminate value of uj. 



4.2. Constraints from recent SNIa observations 

In the recent years the confidence in type la supernovae as standard candles has been steadily growing. Actually 
it was just the SNIa observations that gave the first strong indication of an accelerating expansion of the universe, 
which can be explained by assuming the existence of some kind of dark energy or nonzero cosmological constant 
(Schr mcTt fc al. 1998|l . Since 1995 two teams of astronomers - the High-Z Supernova Search Team and the Supernova 
Cosmology Project - have been discovering type la supernovae at high redshifts. First results of both teams were 
published by Schmidt & al. (1998) and Perlmutter & al. (1999). Recently the High-Z SN Search Team reported 
discovery of 8 new supernovae in the redshift interval 0.3 < z < 1.2 and they compiled data on 230 previously 
discovered type la supernovae ( |Tonry et al. 2001) ). Later Barris & al. (2004) announced the discovery of twenty-three 
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Fig. 10. Behaviour of the CMB power spectrum for the double exponential potential model with Ho ~ 0.97 and u = 0.1 
compared with the data of Boomerang (empty boxes) and WMAP (filled boxes) . 

high-redshift supernovae spanning the range of z = 0.34 — 1.03, including 15 SNIa at z > 0.7 . 

Recently Riess & al. (2004) announced the discovery of 16 type la supernovae with the Hubble Space Telescope. This 
new sample includes 6 of the 7 most distant (z > 1.25) type la supernovae. They determined the luminosity distance 
to these supernovae and to 170 previously reported ones using the same set of algorithms, obtaining in this way a 
uniform gold sample of type la supernovae containing 157 objects. The purpose of this section is to test our scalar field 
quintessence models by using the best SNIa dataset presently available. As a starting point we consider the gold sample 
compiled in Riess & al. (2001). To constrain our models we compare through a x 2 analysis the redshift dependence of 
the observational estimates of the distance modulus, \x — ra — M, to their theoretical values. The distance modulus is 
defined by 

m-M = 51ogL> L (z) + 51og(-^-) + 25, (62) 

-no 

where m is the appropriately corrected apparent magnitude including reddening, K correction etc., M is the cor- 
responding absolute magnitude, and Dl is the luminosity distance in Mpc. For a general flat and homogeneous 
cosmological model the luminosity distance can be obtained through an integral of the Hubble function H , as 

D L (z) = (l + z ) J* J— d (, (63) 



4.2.1. The single exponential potential 



For the single exponential potential model the luminosity distance can be analytically calculated from Eq. Ij63(l , using 
the Hubble function given in Eq. (JHJ, and the z(t) relation as given by Eq. (|15|) . The luminosity distance can be 
represented in the following way: 



D L (t) = (l + z(f)) 



I r[§] 



3( 



2t 



l + t- 



:)* 2*1 



ills 
6' 3' 6' 



(64) 



Inverting the relation z(t), we can construct D^(z), and evaluate the distance modulus according to Eq. (|62|) . 
Performing an x 2 analysis with the gold dataset of Riess & al. (2004) we obtain xled = 1-15 for 157 points, and 
as best fit for H the value H — 0.98lo.o5' which corresponds to f^o = 0. 25^^08' anc ^ ^ = 0-65-o'cm- These values 
are consistent with the WMAP data analyzed above. In Fig. ^] we compare the best fit curve with the observational 
dataset. 

For the double exponential potential model we can perform the same analysis, using Eqs. I|20l) and (|26() . With 
our choice of ui = 0.1 also in this case the z(t) relation can be inverted. Again, through Eq. I|63|l . which in this 
case can be integrated only numerically, we construct the distance modulus and perform the % 2 analysis on the gold 
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Fig. 11. Observational data of the gold sample of SNIa (Riess et al. 2004) fitted to our model with the single exponential 
potential. The solid curve is the best fit curve with Ho = 0.981q'q5. It corresponds to Qmo = 0.25^Q'Qg, and h = 0.65^q'q4. The 
inner rectangle zooms in on the high redshift SNIa. 




Fig. 12. Observational data of the gold sample of SNIa (Riess et al. 2004 ) fitted to the double exponential potential model. 
The solid curve is the best fit curve with H = 0.99±° ™. lt corresponds to f2 M0 = 0.25to° tt ° l, and h = 0.65±£^- The last 
term in the error on Qmo quantifies the effect of the frequency parameter uj, when it changes from uj — 0.1 to uj = 0.9. The 
inner rectangle zooms on the high redshift SNIa. 



data set. We obtain Xred = 1-16 for 157 data points, and the best fit value is Hq = 0.99^o'o7> which corresponds to 
^mo = 0-25^o'o6^o 02: an d h = 0.65loo4- The last set of errors in Omo quantifies the effect of the parameter u>, when 
it changes from uj = 0.1 to uj = 0.9. We noted that for this potential the role of the high redshift supernovae is quite 
important, and they change the value of Xred from 1.1, if we consider the supernovae at z < 1, to Xred = 1-17, once 
we use the whole data set. This circumstance confirms the necessity to increase the statistics of data at high redshifts 
to discriminate among different models. In Fig. 1121 we compare the best fit curve with the observational data. 
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4.3. Constraint from galaxies redshift surveys 

Once we know how the growth index / changes with redshift and how it depends on CIm we can use the available 
observational data to estimate the present value of SImo- The 2dFGRS team has recently collected positions and 
redshifts of about 220,000 galaxies and presented a detailed analysis of the two-point correlation function. They 

/ 

measured the redshift distortion parameter 3 = — , where b is the bias parameter describing the difference in the 

o 

distribution of galaxies and mass, and obtained that /3| 2 ^o.i5 = 0.49±0.09 and b = 1.04±0.11. From the observationally 
determined [3 and b it is now straightforward to get the value of the growth index at z = 0.15 corresponding to the 
effective depth of the survey. Verde & al. (2001) used the bispectrum of 2dFGRS galaxies, and Lahav & al. (2002) 
combined the 2dFGRS data with CMB data, and they obtained 

b verde = 1.04 ±0.11, (65) 
blahav = 1.19 ±0.09. (66) 
Using these two values for b we calculated the value of the growth index / at z = 0.15, we get respectively 

h = 0.51 ±0.1, (67) 
f 2 = 0.58 ±0.11. (68) 

4.3.1. The single exponential potential 

Using Eq.(15) we express time t through the redshift z in Eq. I|34|) and then substituting z = 0.15 and the two values 
of fi and fi we calculate Hq. Substituting the thus obtained Hq into Eq.(9) and setting t = 1 we get 



Qmoi = 0.26 ±0.1, 
Q M02 = 0.32 ±0.13, 



(69) 
(70) 



which can be weighted to give the final value ^mo = 0.27 ±0.09. Substituting the previously obtained relation between 
time t and the redshift z and Hq into Eq.(9) we calculate £Im(z = 0.15) and we get 



fi M i = 0.32 ±0.1, 
n M2 = 0.39 ±0.13, 



(71) 
(72) 



which gives the final value £Im(z = 0.15) = 0.35 ± 0.09. It turns out that this value is also fully compatible with the 
independent estimates derived, for instance, from the first data release of the SDSS (SDSS colla boration 2003~j) . 



4.3.2. The double exponential potential 

Repeating the same procedure in the case of models with double exponential potential (with u> — 0.1), we get 

Omoi = 0.26 ±0.1, (73) 

f7 M0 2 = 0.31 ±0.11, (74) 

which can be weighted to give the final value f^Mo = 0-28 ± 0.07 ± 0.04, where the last set of errors in Omo quantifies 
the effect of ui when it changes from to = 0.1 to uj = 0.9. As in the previous case, using Eq. I|27|l we express time t 
through the redshift z and using Eq. (21) we calculate CIm{z — 0.15); we get 



Q M1 = 0.32 ±0.1, 
Q M2 = 0.38 ±0.11, 

which gives the final value n M (z = 0.15) = 0.35 ± 0.07 ± 0.04. 



(75) 
(76) 



4.4. Comparison with the standard AC DM model 

In the standard ACDM model we can write the growth index as a function of the redshift z and f^MOi obtaining: 



/a = 



{l + zf 2*1 


1 1 11 ( l-O.Mo\ 

3' x ' 6 ' ^ n M0 J 




2F1 


1 1 11 ( 1-fiMO ^ 
[a' 4 "' 6 ' V°Mo(l+z) a y 







(77) 



2-F1 


1 1 11 ( 1-fiMO ^ 

3' x ' 6 ' \^n M o ) 






{l + zf 2 F 1 


1 1 11 ( 1-^MO^ 
3' ±5 6 ' ^ f2 M o J 





6(1-Omo) 2 ^i[|,2,¥,-(^Jm^) 
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Fig. 13. The growth index / in several cosmological models: the solid line corresponds to the ACDM model with Omo = 0.25. 
The dashed curves correspond to the two exponential potentials: the thin line is for the model with single exponential potential 
and the thick one is for the double exponential potential model. 



Applying the same procedure as in the case of the two quintessence models we get 

fl M01 = 0.22 ±0.11, (78) 
n M 02 = 0.28 ±0.13, (79) 

and a final estimate of ^mo — 0.25 ±0.09, which is compatible with the ^mo obtained from independent measurements 
(see for instance IPerlmutter fc al. 199911 and our estimates. Therefore using only the value of Qmo it is not possible 
to discriminate between the models with the standard cosmological constant and the quintessence models considered 
here. However, as is seen in Fig. 1131 independent measurements from large redshift surveys at different depths can 
disentangle this degeneracy. 



5. The dark energy and the Hubble age 

It is known that the age estimates of the universe can in principle independently constrain the equation of state of 
the dark energy, since the age of the universe, as well as the lookback time - z relation, is a strongly varying function 
of w QKrauss 20041 iFerreras fcTi ilk 2003). In the approximation of constant w, the age of a flat universe is defined by 
the formula: 

H t = — ; r-r, (80) 

where Qx is the fraction of the closure density in material with the dark energy equation of state. The lookback time 
is the difference between the present age of the universe and its age at redshift z, and it is given by: 

t(z)=t H f — ; r-r, (81) 

where tn is the Hubble time. It is worth to note that since the lookback time - z relation does not depend on the 
actual value of to, it furnishes an independent and interesting cosmological test trough age measurements, especially 
when it is applied to old objects at high rcdshifts (Alc aniz fc al. 2003}) . In principle it is possible that a model fits well 
the lookback time data and at the same time gives a wrong value for to- For varying w, as in the case of our models, 
the Eq. H80|) can be rewritten in a more general form as 

Jo aF ia) 

where F(a) = Sg^, while the lookback time relation becomes 



(83) 
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On 

Fig. 14. The curve t(w, Q M ) = 1. 

The purpose of this section is to evaluate the Eqs. (|82[1 and H83[) for our quintessence models in order to constrain the 
value of the parameter w. Before going further we note that this analysis is at the same time simple and particularly 
interesting in the context of our parametrization, which is based on the choice of the present age of the universe as a 
unit of time. 



5.1. The single exponential model 

As a starting point we note that since the age of the universe has been set equal to unity (to = 1), the dimensionless 
quantity Hgto which appears in Eq. (|82() is simply Hq. This means that Hq is the only parameter that enters in the 
observational quantities. From our previous analysis we obtained H = 0.97, which agrees with the limit of H t < 1.1 
obtained in Krauss (2004). Let us start with the lookback time - z relation, which has the form 



1 



2-3H, 



2- 



9 Hi (1 + zf + 8 z (z 2 + 3z + 3) (2 - 3 H ) 



(1 



Using the relation 

(A-3H )((3H -2)t 2 



3#o 



[2(3iJ -2)i 2 +4-3H ] 2 
we obtain t = t(z, £Im) as 

V2 



t{z,Q 



M ) 



(1 + z)y/y/\+m M -(l + 2Sl M ) 



(84) 



(85) 



(86) 



{Vl + 8fi A/ (1 + zf - (1 + AVL M ) (1 + 2z) + z 2 (1 - 40 M ) 
-a/2 (l + z) Ul + 8il M (1 + (1 + 4fi M ) (3z + 3z 2 + z 3 )) 

1 

- (1 + 4fi M ) - z (3 + 3z + z 2 ) (1 + 80 M (1 + ^m))] * } 2 • 
If we impose, according to our assumptions, that i(0, fijvf) = 1 we obtain SImq = 1 — 



(37?o - 2) (3H + 4) 



9H* 



which, 



for Hq = 0.97 ± 0.08 gives Hmo = 0.26 ± 0.05. In a similar way we can construct the surface t(w, fW), eliminating z 
between Eq.JHH) and (13). The curve t(w, fijw) = 1 determines the physically acceptable values of parameters in the 
plane w — flM, as shown in Fig ll4l 

In Fig ll5l we show the dependence of the dimensionless quantity -f/o^o on the value of the parameter w. 



5.2. The double exponential potential 

For the model with double exponential potential the procedure outlined for the single exponential potential becomes 
more complicated from the computational point of view. Actually the z — t relation in Ea. l|26() is rather involved and 
cannot be exactly inverted. However, since t € (0, 1) and for small values of u>, it can be simplified to the invertible 
form of Eq. Ij27(l . by means of a series expansion. In Fig 1 161 we compare the exact result with the approximate one. 
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Fig. 15. Dependence of the dimensionless quantity Hoto on w for the single exponential model. 




Fig. 16. Behaviour of the exact and approximate t — z relation for the double exponential model. The solid line is the exact 
function and the dashed one the approximate function. 




Fig. 17. Dependence of the dimensionless quantity Hoto on w, for the model with double exponential potential. 



Once t(z) is known we can apply the same procedure as in the single exponential case: it is possible to construct 
the surface t(w, fijvf), eliminating z between the w(t), and CIm- The curve t(w, Ojf) = 1 determines the physically 
possible values of parameters in the plane w — £Im . 

In FigEl we show the dependence of the dimensionless quantity Hoto on the value of the parameter w. 

6. Conclusions 

We tested the viability of two classes of cosmological models with scalar field models of quintessence with exponential 
potentials. To compare predictions of our models with observational data we used three different types of observations: 
namely, the power spectrum of the CMB temperature anisotropy and in particular positions and heights of the peaks, 
the high-redshift SNIa data compiled in Riess & al. (2004) (gold data set), and the large survey of galaxies by the 
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2dFGRS team (Hawkin s fc al. 2003| ) and in particular their estimate of the average density of dark matter. We showed 
that predictions of our models 

Table 1. The basic cosmological parameters derived from our models are compared with those presented by the 
WMAP team. The position of the third peak Z3 has been determined by the Boomerang team. The last set of errors 
in the case of the double exponential potential quantifies the effect of the parameter ui, when it changes from oj = 0.1 
to oj = 0.9. 







(A exp ( \ ypy) - B exp (- § y/\<p)) " 


WMAP 




0.13 ±0.03 


0.13 ±0.03 ±0.04 


0.135 ±0.008 




0.75 ±0.08 


0.76 ± 0.07 ± 0.04 


0.73 ±0.04 


w 


-0.76 ±0.05 


-0.74 ±0.05 ±0.03 


w < -0.78 


r 


1Q+ UU7 


n , Q +u.u7+u.ua 

u - 1;, -0.05-0.02 


0.19 ±0.02 


R 


-1 71 +0.U9 
L - 1 L -0.09 


1 71 +U.U9 
L - 1 L -0.09 


1.71 ±0.137 


h 




224+g 1 


220.1+- 


h 


546+^ 


543±iS 




h 


8301S 


825^34 


851±|i 



are fully compatible with the recent observational data. However the models considered here, despite being fully 
compatible with the present day observational data have very different behavior on long time scales. The model with 
the single exponential potential will eternally accelerate, and though for sufficiently large t its evolution is dominated by 
dark energy with VL V — + 1 and Qm — > the scale factor increases only as a power law a(t) ~ t 4 ^ 3 and not exponentially. 
In this model there is a particle horizon. The model with the double exponential potential has a different asymptotic 
(t — > 00) behavior. First of all, though a(t) — > 00 as t — > 00, Qm does not decrease to zero, so for large time matter 
and dark energy coexist. For sufficiently large t the scale factor grows like in a FRW matter-dominated flat model 
a(t) ~ t 2 / 3 . In this model, the particle horizon does not form and asymptotically ( for large t) the expansion of 
the universe is decelerating. It also turns out that the available observational data cannot discriminate between the 
quintessence models considered here and the ACDM model. New data coming from high redshifts observations could 
remove this degeneracy. 
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